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Abstract 

This  survey  is  a  summary  of  available  information  on  the  corpuscular  radia¬ 
tion  environment  in  the  vicinity  of  the  earth  and  in  the  earth's  atmosphere;  radia¬ 
tion  trapped  in  the  magnetosphere  (Van  Allen  belts),  galactic  cosmic  rays,  and 
high~energy  solar  particles  (solar  cosmic  rays)  are  described.  The  topics 
presented  are;  the  motions  of  charged  particles  in  the  geomagnetic  field;  the 
spatial  and  temporal  distributions,  intensities,  and  energies  of  electrons  and 
protons  trapped  in  the  earth's  magnetic  field;  the  composition,  intensity,  energy 
spectrum  and  variations  of  galactic  cosmic  rays;  and  the  composition,  intensities, 
anisotropy,  and  variations  of  the  high-energy  charged-particles  emitted  sporadi¬ 
cally  by  the  sun  in  association  with  solar  flares  and  other  solar  disturbances. 
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This  report  is  a  revision  of  Sections  17.  2  through  17.  5  of  Chapter  17 
"Corpuscular  Radiation"  of  the  Handbook  of  Geophysics  and  Space  Environments’*' 
changes  in  Section  17.  1  are  minor.  (Numbers  of  Sections  are  the  same  as  those 
in  the  original  Handbook  ,-o  that  the  cross-referencing  system  in  other  chapters 
remains  valid. )  This  survey  represents  the  state  of  the  art  in  June,  1968,  when 
the  manuscript  was  submitted. 
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Scientific  Editor 

Handbook  of  Geophysics  and  Space  Environments 


*  Published  by  the  Air  Force  Cambridge  Research  laboratories  and  by  the 
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Corpuscular  Radiation  :  A  Revision  of 
Chapter  17,  Handbook  of  Geophysics 
and  Space  Environments 


This  chapter  outlines  current  knowledge  of  corpuscular  radiation  in  the  vicinity 
of  earth;  this  includes  the  charged  particles  trapped  in  the  geomagnetic  field,  cos¬ 
mic  radiation,  and  high-energy  soiar  particles  associated  with  solar  flares.  The 
sources,  acceleration  mechanisms,  diffusion  through  space,  and  solar  modulation 
processes  are  not  well  understood.  The  prediction  of  the  occurrence  of  high-energy 
solar  particle  emission  and  of  the  characteristics  of  corpuscular  radiation  emitted 
during  solar  flare  events  is  still  aty  unsolved  problem. 

17.1  B  VSlt:  RKUTIONS 

17.1.1  esrUuwe  sad  tMlsW— » 

Corpuscular  radiation  is  usually  characterised  by  the  energy  of  the  particle  in 
electron  volts.  (The  international  symbol  tor  109  eV,  now  adopted  by  the  U.S. 
Bureau  of  Standards,  is  GeV,  rather  than  BeV. )  The  total  energy  is  the  sum  of 
the  kinetic  energy,  E^,  plus  rest-mass  energy, 

Ef  •  Ejj  ♦  n»0C2  ■  me2,  (17-1) 


(Received  for  publication  2  December  19fi8) 


where  m  is  the  rest  mass,  m  the  relativistic  mass,  and  c  the  speed  of  light  in 
o  n  . 

vacuo.  For  the  proton,  mQc"  is  938  MeV;  for  the  electron,  mQc  is  0.  511  MeV. 

It  is  customary  in  the  literature  to  give  values  of  the  energy  per  nucleon  for 

heavier  particles. 

The  total  energy  is  related  to  the  particle  momentum,  p,  by 


Et  *  (p  c 


2  2  .  ...  2  .4,1/2 


+  m  “c  ) 
o 


(17-2) 


In  the  absence  of  electric  fields  and  of  time-varying  magnetic  fields,  and  p  are 
constants  of  the  particle's  motion.  In  cosmic  ray  physics  the  momentum  is  often 
given  in  units  of  eV/c. 

Observations  of  the  intensity  of  corpuscular  radiation  are  reported  in  various 
units  depending  upon  the  detection  method.  The  unidirectional  differential  intensity, 
J(E),  is  the  flux  (number  per  unit  time)  of  particles  of  a  given  energy  per  unit 
energy  interval  in  a  unit  solid  angle  a.  >ut  thr  direction  of  observation,  incident  on 
a  unit  area  perpendicular  to  the  direction  of  observation;  the  units  are  usually  par- 
tides  cm  s  sr  MeV  .  Unidirectional  integral  intensity,  J  (>  E),  is  the  in¬ 
tensity  of  particles  with  energy  greater  than  a  threshold  energy,  E, 


J  (>E)  = 


J(E)dE. 


(17-3) 


Omnidirectional  intensities  are  J(E),  or  J(>E),  integrated  over  4tt  solid  angle. 

The  diffei  mtial  energy  spectrum  is  J(E)  plotted  against  E.  The  integral 
energy  spectrum  is  J  (>E)  plotted  against  E.  J(E)  and  J(>E)  sometimes  can  be  ex¬ 
pressed  as  a  power  law; 

J  (E)  *  JqE"7  ,  (17-4) 

J  (>E)  «  JQ(y  -  1)'1E'Y+1  .  (17-5) 


where  y  >  1.  It  is  also  convenient  sometimes  to  use  an  exponential  law; 


J(E)  ■  JQ  exp  (-r/EQ), 


(17-6) 


J(>E)  ■  JQ  EQ  exp  (-E/E0) 


(17-7) 
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where  Eq  is  defined  as  the  average  particle  energy  if  the  spectrum  is  of  this 
simple  exponential  form. 

Observations  of  counts  per  unit  time  or  of  ionization  rates  arc  often  difficult 
to  interpret  and  convert  to  J  (E)  or  J  (>E).  V  .'hen  considerations  of  the  effects  of 
corpuscular  radiation  on  personnel  are  of  primary  concern,  intensities  or  counting 
rates  must  be  converted  to  dose  or  dose  rates.  The  rad  is  the  unit  of  absorbed 
dose;  one  rad  is  JO0  erg  absorbed  per  gram  of  absorbing  material.  Conversion  of 
observed  particle  fluxes  to  absorbed  dose  rate,  however,  is  not  straightforward  and 
is  sometimes  impossible  for  lack  of  necessary  details;  the  conversion  depends  in 
a  complex  way  on  the  energies  and  kinds  of  particles  and  on  geometrical  configura¬ 
tions  of  the  absorbers  and  direction  of  the  incident  beam,  as  well  as  on  the  different 
absorbing  properties  of  materials. 

The  magnetic  rigidity.  P,  of  a  particle  is  a  measure  of  its  resistance  to  a 
magnetic  force  that  deflects  the  particle  from  a  straight-line  trajectory.  The 
rigidity  is  defined  as 

P  =  pc/q  .  (17-8) 

where  q  is  the  charge  of  the  particle.  If  pc  is  in  electron  volts,  then  q  is  the  num¬ 
ber  of  electronic  charge  units  and  P  is  in  volts. 

Figure  17-1  shows  the  relation  between  the  rigidity  of  protons  and  alpha  par¬ 
ticles  and  their  kinetic  energy  per  nucleon;  in  using  this  figure  for  conversion  of 
units,  note  that  the  energy  scale  must  be  multiplied  by  4  (the  number  of  nucleons) 
to  obtain  the  total  kinetic  energy  of  an  alpha  particle. 

On  many  occasions  the  magnetic  rigidity  is  used  in  place  of  the  kinetic  energy 
in  describing  the  flux  spectrum.  In  particular,  in  lieu  of  Eq.  17 -t>,  a  frequently 
used  expression  is: 

J(P)  -  J0  exp  (-P/PQ)  (17-9) 

The  geomagnetic  field  acts  on  cosmic  rays  as  a  magnetic  analyzer,  removing 
low  rigidity  charged  particles  from  the  flux  incident  at  the  top  of  the  earth's  atmos¬ 
phere;  particles  of  progressively  lower  rigidity  are  detected  as  the  distance  north 
or  south  of  the  geomagnetic  equator  increases.  The  cutoff  or  threshold  rigidity  of 
cosmic  rays  is  the  minimum  rigidity  that  permits  a  charged  particle  to  arrive  from 
a  specified  direction  at  a  given  latitude  and  longitude;  particles  of  lower  rigidity 
are  not  observed  at  the  specified  location  and  direction  because  of  this  geomagnetic 
cutoff. 


ENERGY  PER  NUCLEON  (GeV) 
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Figure  17-1.  Conversion  From  Magnetic  Rigidity  to  Kinetic  Energy  per 
Nucleon  for  Protons  and  Alpha  Particles.  Because  nuclei  heavier  than 
He3  have  about  the  same  charge  to  mass  ratio  as  alpha  particles  (0.  5). 
the  alpha -particle  curve  also  gives  approximate  values  for  the  heavier 
nuclei 


At  a  given  location  on  the  surface  of  the  earth,  only  those  cosmic  ray  particles 
that  arrive  from  certain  specific  directions  with  respect  to  the  station  can  be 
detected.  These  allowed  directions  form  the  asymptotic  cone  of  acceptance;  for 
asymptotic  cones  of  acceptance  of  various  cosmic-ray  stations  see  McCracken, 
et  al  (1965). 
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17.1.2  Trapped  Particle  Motion  in  a  Magnetic  Field 
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Figure  17*2.  Diagram  of  the  Simplified 
Motion  of  a  Charged  Particle  Trapped  in  a 
Magnetic  Dipole  field;  Azimuthal  Drift  is 
Ignored.  V  is  the  particle  velocity,  if  the 
magnetic  vector  of  the  dipole  field,  and  a  is 
the  pitch  angle.  (After  Singer  and  Lenchek, 
1962. ) 


Figure  17-2  illustrates  the 
motion  of  a  charged  particle  trapped 
in  a  magnetic  dipole  field.  The  par¬ 
ticle  is  confined  to  the  surface  of  a 
tube  of  magnetic  lines  of  force,  mov¬ 
ing  in  a  quasi-helical  path  with  a  ro¬ 
tational  period  t^,  and  oscillating 
back  and  forth  between  mirror 
points  with  a  period  The  pitch 
angle,  a,  is  the  angle  between  the 
momentum  vector,  p,  and  the  mag¬ 
netic  vector,  B;  at  the  mirror 
point,  a  is  90°.  The  particle  also 
drifts  slowly  around  the  earth  with 
an  azimuthal  drift  period  r3;  nega¬ 
tive  particles  drift  to  the  east, 

positive  particles  to  the  west.  In  the  geomagnetic  field  the  orders  of  magnitude  of 

*6  “3 

the  periods  are:  for  electrons,  10  s;  for  protons,  10  s;  for  electrons 

and  protons,  1  s;  for  electrons  and  protons,  103  s  . 

The  period  r ^  is  the  reciprocal  of  the  gyrofrequency  (cyclotron  frequency). 

Tj  is  27rm/qB  or  2v  ET/c2qB  in  System e  International  (MKSA);  in  Gaussian  (mixed 

is  p(sina)/qB;  in  Gaussian 

Because  geomagnetic  lines  of  force  converge  toward  the 

poles,  the  particle  penetrates  regions  of  increasing  magnetic  flux  density  as  it 

leaves  the  equatorial  region.  It  is  repelled  by  a  force  that  depends  upon  r  and 

6 

gradient  B,  so  that  the  component  of  momentum  parallel  to  the  magnetic  line  of 
force,  p  cos  a,  decreases  and  becomes  zero  (a  ■  90°)  at  the  mirror  point.  The  par¬ 
ticle  then  spirals  back  to  regions  of  lower  magnetic  flux  density.  In  the  absem  e  of 
an  electric  field,  the  magnitude  of  the  momentum,  |  p  j ,  is  constant. 

If  the  spatial  and  temporal  variation  of  the  magnetic  field  is  not  too  rapid,  so 
that  Tg  | grad  B  |«  B  and  TjdB/dt«B,  the  path  can  be  considered  a  circular  mo¬ 
tion  about  a  center  that  moves  along  the  line  of  magnetic  force  perpendicular  to  the 

circle  and  located  at  its  center;  this  center  is  the  guiding  center.  The  rotating 

2  2  o 

charge  generates  a  magnetic  moment  parallel  to  B  and  of  magnitude  p  c  (sin"a)/2ETB, 
The  repelling  force  is  the  product  of  this  magnetic  moment  and  dB/ds,  where  ds  is 
the  element  of  length  along  the  guiding  center  path.  Under  the  conditions  of  slow 
variation  for  which  the  guiding  center  approximation  is  valid,  the  magnetic  moment 
is  constant;  it  is  called  the  first  adiabatic  invariant.  For  a  constant  magnetic  moment 

O 

(constant  magnetic  flux  through  the  area  nr^  ), 


CGS)  units,  Tj  is  27rET/cqB.  The  radius  of  gyration,  r  , 


units,  rg  is  pc(sina)/qB. 
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Bm=B/sin2a,  (17-10) 

where  B  is  the  magnetic  field  at  a  position  where  the  pitch  angle  is  or,  and  is 
the  field  at  the  mirror  point.  Bm  depends  only  on  the  initial  conditions  of  the 
particle  direction  and  the  magnetic  field;  it  is  independent  of  particle  charge,  mass, 
and  energy. 

The  longitudinal  invariant,  or  second  adiabatic  invariant,  is 


J 


a 


s2 


(p  cos  a)  ds  , 


(17-11) 


where  ds  is  the  element  of  path  length  of  the  guiding  center  from  one  mirror  point, 
si,  to  its  conjugate,  s2.  From  Eq.  (17-10)  and  (17-11) 

s2 

1  =  Ja/p  -  f  [  1  -  (B/Bm)]  1/2  ds  .  (17-12) 

si 

I  is  the  invariant  used  in  establishing  the  magnetic  shell  parameter  L  (Sec.  17.  2.  1); 
if  the  guiding  center  approximation  is  valid  and  if  r 2  j  dB/dt  |  «E,  I  is  conserved. 

The  bounce  period,  r 2,  is  the  time  of  travel  from  one  mirror  point  to  the 
other  and  back; 


T 


2 


■  2 


si 


E^ds 

~2 

c  p  cos  a 


(17-13) 


If  changes  in  the  magnetic  field  are  small  within  one  azimuthal  drift  period, 
then  there  exists  a  third  adiabatic  inva  ”iant  called  the  flux  invariant.  W  ithin  the 
guiding-center  approximation,  the  flux  invariant  is  equal  to  the  magnetic  flux  en¬ 
closed  by  the  guiding-center  trajectory  during  a  single  drift  around  the  earth.  For 
example,  if  the  geomagnetic  field  is  slowly  compressed  (or  expanded)  the  particle 
will  also  move  inward  (or  outward)  in  such  a  way  that  the  magnetic  flux  enclosed 
by  the  guiding-center  trajectory  will  remain  constant. 
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17.2  TRAPPED  RADIATION 


For  convenience  of  discussion  and  consistency  with  definitions  in  the  litera¬ 
ture,  the  trapped  electrons  and  protons  are  classified  according  to  their  location 
as  follows: 


(a)  The  inner  zone  (also  called 


the  inner  Van  Allen  belt)  where 
large  temporal  variations  of  the  mag¬ 
netic  field  do  not  occur  and  particle 
lifetimes  are  relatively  long;  and 
(b)  The  outer  zone  (also  called 
the  outer  Van  Allen  belt)  where  sig¬ 
nificant  temporal  variations  of  the  geo- 
matnetic  field  can  occur,  particle  life¬ 
times  are  relatively  short,  and  solar 
produced  particles  may  be  trapped  (at 
least  temporarily).  A  different  classi¬ 
fication,  according  to  location  in  the 
magnetosphere,  is  becoming  increas¬ 
ingly  important  as  trapped  radiation  is 
considered  a  part  of  magnetospheric 
phenomena.  Figure  17-3  illustrates 
the  regions  of  stable  trapping,  pseudo 
trapping,  and  no  trapping  in  the  mag¬ 
netosphere.  Most  of  the  trapped 
radiation  is  contained  in  the  region  of 
stable  trapping,  but  during  disturbed 
geomagnetic  conditions  the  pseudo 
trapping  region  may  contain  trapped 
particles  which  will  be  lost  in  the 
earth's  atmosphere  at  high  latitudes. 
The  earliest  measurements  of 


60"N  90°N  60°N 


Figure  17-3.  Cross  Section  of  the  Noon- 
Midnight  Meridian  Plane  of  a  Model 
Magnetosphere,  Showing  Computed  Limits 
of  Stable  Trapping  and  Pseudo  Trapping 
Regions.  Particles  mirroring  inside 
pseudo-trapping  regions  are  unable  to 
complete  a  180®  drift  around  the  earth; 
those  injected  into  the  left-hand  side  will 
be  lost  into  the  tail  whereas  those  injected 
into  the  right-hand  side  will  abandon  the 
magnetosphere  on  the  day  side  through 
the  boundary.  (After  Roederer,  1967). 


trapped  particles  were  obtained  with 

relatively  simple  detectors  that  could  not  distinguish  uniquely  between  various  par¬ 
ticle  species  and  energies;  hence  the  interpretation  of  these  early  data  is  somewhat 
ambiguous.  Since  late  1962,  the  instruments  flown  have  been  sufficiently  sophisticated 
to  enable  unique  Interpretation  of  particle  species  and  energy  in  all  regions  of  the 


magnetosphere. 
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17.2.1  Spill  ml  l>tirtimi«lrrn  ('huriu  tcrizin^  Trapped  Purliclos 

In  order  to  determine  and  describe  the  trapping  regions,  a  two-parameter 
convention  has  evolved,  which  is  based  on  the  motion  of  the  guiding  center  along  a 
magnetic  field  line.  In  developing  these  parameters,  the  rapid  particle-motion 
perpendicular  to  the  field  is  ignored.  If  the  slow  drift  around  the  earth  is  also 
ignored  the  guiding  center  moves  along  a  single  line  of  force  between  the  mirror 
points.  (The  drift  can  be  treated  as  a  perturbation  of  this  guiding-center  motion. ) 
Figure  17-2  shows  the  guiding-center  path,  the  lines  of  magnetic  force,  and  the 
reflection  from  one  hemisphere  to  the  other. 

For  any  given  azimuthal  angle,  two  parameters  suffice  to  describe  this  guiding- 
center  motion.  Two  such  parameters  (by  no  means  the  only  possible  ones)  are 
B^,  the  scalar  magnetic  field  at  the  point  where  the  particle  is  reflected,  and  Rq, 
the  maximum  excursion  of  the  guiding-center  trajectory  (i.  e.  the  maximum  dis¬ 
tance  of  the  magnetic  field  line  from  the  center  of  the  geomagnetic  field. )  Because 
B  is  an  adiabatic  invariant  of  the  particle  motion,  it  is  a  useful  parameter  with 
an  obvious  physical  meaning.  However,  Rc  is  not  an  invariant.  Due  to  the 
azimuthal  asymmetry  of  the  geomagnetic  field,  for  a  given  particle  can  vary  as 
much  as  10%  as  the  particle  drifts  around  the  earth.  The  variation  depends  not 
only  on  the  initial  value  of  Rq  and  the  azimuthal  angle,  but  also  on  B  . 

The  second  adiabatic  invariant,  the  longitudinal  invariant  I  (Eq.  17-12),  is 
available,  but  this  does  not  have  a  clear  physical  or  geometrical  meaning.  In  a 

Q 

pure  dipole  field  with  moment  M.  Rq  is  determined  by  Bm  and  I;  Ro  Bm/M  * 
f(l3Bm/M).  In  the  geomagnetic  field  this  relation  does  not  hold.  Mcllwain  [  1081) 
replaced  Rq  by  L  and  defined  this  new  parameter  by  the  relation 

L3Bm/M  -  f(I3Bm/M)  .  (17-  H) 


where  the  units  of  L  are  earth  radii,  (Rf.),  and  M  »  0.  311853  gauss-Rj,.3. 

For  particles  traveling  along  a  given  field  line,  the  L  value  depends  only 
slightly  on  where  the  particles  mirror.  Because  the  geomagnetic  field  is  nearly  a 
dipole,  L  %  Rq,  and  L  has  an  approximate  geometrical  meaning.  L  can  be  cal¬ 
culated  for  any  part  of  the  geomagnetic  field  of  interest;  therefore,  L  and  B  are 
used  extensively  to  describe  the  trapping  regions. 

It  is  implicit  in  the  definition  of  L  that  the  value  of  I,  at  any  given  point  in 
space  is  an  adiabatic  invariant  of  the  motion  of  those  particles  which  mirror  at  the 
given  point  |  Stone,  1663] .  All  particles  passing  through  the  point  do  not  in  general 
have  that  particular  value  of  L,  but  they  do  have  L  values  within  a  few  percent  of 
it.  The  locus  of  all  points  with  a  given  1.  value  is  called  an  L-shell.  All  particles 
with  a  particular  I.  must  mirror  on  the  corresponding  L-shell,  but  their  motion  may 
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otherwise  deviate  from  the  siiell  by  distances  equal  to  a  few  percent  of  L.  From  a 
practical  standpoint,  the  accuracy  of  measurements  has  not  progressed  to  a  point 
where  these  deviations  from  an  L-shell  are  significant,  hut  ignoring  them  lias  led 
to  confusion  as  to  the  significance  of  l,.  As  more  precise  data  become  available, 
these  details  assume  greater  importance.  Ftgury  17-4  shows  contours  of  constant 
I.  and  11  at  zero  geographic  longitude.  Values  of  the  magnitude  of  11  at  various 
altitudes  and  locations  are  given  in  Chapter  11;  these  values,  however,  should  not 
he  used  for  estimating  fluxes  below  1000 -km  altitude  in  the  South  Atlantic  region. 

Care  must  he  exercised  in  comparing  data  on  trapped  particles  when  different 
models  of  the  geomagnetic  field  were  used  to  calculate  the  H-I.  coordinates.  There 
arc  significant  differences  between  various  geomagnetic  models,  particularly  at 
low  altitudes  over  regional  anomalies,  that  can  result  in  apparent  differences  of 
trapped  particle  fluxes  in  11-1.  space;  for  a  discussion  of  this  problem,  see 
l.indstrom  and  Heckman  (1997). 

The  pitch  angle  at  the  geomagnetic  equator  is  used  sometimes  as  a  more  con¬ 
venient  parameter  than  11  for  interpreting  results.  In  practice,  the  value  of  I,  is 
often  calculated  at  the  equator.  The  equatorial  value  does  not  correspond  exactly 
to  l.  at  the  mirror  point,  hut  is  accurate  to  a  few  percent. 

A  coordinate  system  u-eful  in  obtaining  an  intuitive  feeling  for  experimental 
results,  as  well  as  in  the  interpretation  of  data,  effectively  transforms  the  real 
geomagnetic  field  into  a  pure  dipole  representation.  This  is  accomplished  by  sub¬ 
stituting  for  11  and  I,  in  the  equations  for  a  pure  dipole  field,  the  11  and  1,  of  the 
real  geomagnetic  field.  The  coordinates  are:  R  -  the  distance  from  the  center 
of  an  effective  dipole  (not  the  centered  dipole  defined  in  Section  1 1.  ‘J.  12.  1);  and 
A  -  the  equivalent  geomagnetic  latitude  (similar  to,  but  not  identically  the  geo¬ 
magnetic  latitude  used  in  Section  17. 3). 

I'he  R,  '  coordinates  are  obtained  for  any  point  in  real  space  hv  calculating 
values  of  11  anil  I,  for  the  main  geomagnetic  field,  and  substituting  these  values  in 
the  following  equations: 


R  + 


(3M“/l.B2)  R  -  4M"/ll”  =  0  . 


(17- if,) 


and 


cos2  A  =  R/l.  , 


(17- If.) 


where  R  and  I,  are  in  earth  radii  and  M  is  0.  31 1  (1  f» :i  gauss- R j,*  .  Kquation  (17-15) 
must  be  solved  numerically,  hut  when  R  and  A  are  specified,  I.  is  obtained  directly 
from  Kq.  (17- 1C),  and  11  from 


11  -  M  (4  -  3  cos2  A)1/2/!!2  . 


(17-17) 
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Figure  17*5  shows  the  relation  between  parameters  Bm  and  L  and  the  coordinates 
R  and  X.  Essentially,  these  coordinates  distort  the  geographic  coordinates  so  that 
the  earth  becomes  irregularly  shaped,  and  the  surface  R  >  1  in  this  coordinate 
system  (referred  to  as  the  invariant  earth)  is  not  the  earth's  geoid.  The  actual 
earth's  surface,  the  invariant  surface,  and  a  1000-km  polar  orbit  for  0°  longitude 
in  the  R.X  system  are  included  in  Figure  17-14  for  comparison. 

17.2.2  Natural  Ujectioo  aad  Loss  of  Particles 

The  mechanism  by  which  protons  and  electrons  are  injected  into  the  trapping 
region  is  not  yet  established.  The  hypothesis  of  a  neutron-albedo  source  explains 
the  presence  of  some  of  the  trapped  high-energy  protons.  High-energy  neutrons 
are  created  in  the  upper  atmosphere  by  cosmic  rays,  and  some  of  the  neutrons 
generated  at  the  top  of  the  atmosphere  will  be  reflected  back  into  space;  these 
reflected  neutrons  are  called  the  neutron  albedo.  A  free  neutron  has  a  mean  life¬ 
time  of  about  12  minutes;  it  decays  into  a  proton  with  almost  the  energy  of  the 
neutron,  an  electron  with  any  energy  up  to  725  keV,  and  an  antineutrino.  If  the 
neutron  decays  in  the  trapping  region  of  the  magnetosphere,  the  high  energy 
proton  (and  the  electron)  can  be  trapped.  This  neutron-albedo  source,  however, 
is  small  and  inadequate  to  supply  the  observed  particle  intensities  (see  Fig.  17-20). 
A  possible  source  of  trapped  particles  is  solar  corpuscular  radiation  that  diffuses 
across  magnetic  field  lines  and  is  then  accelerated.  Lack  of  sufficiently  precise 
observations  at  present  has  prevented  the  checking  of  several  theoretical  models 
proposed  for  this  type  of  source. 

Two  important  mechanisms  by  which  high-energy  trapped  protons  escape  are 
interactions  with  the  atmosphere  and  the  failure  of  one  or  more  of  the  adiabatic 
invariants.  Atmospheric  interaction  losses  generally  occur  In  the  region  in  which 
a  particle  mirrors.  Collision  with  atmospheric  particles  decreases  the  trapped 
particle's  energy  by  ionu  .ion  loss.  Small-angle  scattering  lowers  the  mirror 
point  of  any  particle  that  is  at  or  very  near  its  mirror  point.  Catastrophic  loss  is 
a  collision  in  which  the  trapped  particle  scatters  directly  into  the  atmosphere  or 
loses  most  of  its  energy.  Loss  in  the  outer  tone  appears  to  be  caused  by  failure 
of  the  trapping  mechanism.  This  failure  can  result  from  the  deviation  of  the  geo¬ 
magnetic  field  from  a  pure  dipole,  or  from  geomagnetic  disturbances  due  to  solar 
effects  (Sec.  11. 4. 1).  The  breakdown  of  the  trapping  mechanism  is  presumed  to 
limit  the  energy  of  a  proton  that  can  be  trapped  for  long  periods  at  high  altitudes. 
Figure  17-6  shows  the  relation  between  maximum  proton  energy  end  the  magnetic 
shell  parameter  as  obtained  from  experimental  data  and  as  calculated  by  considering 
breakdown  of  the  trapping  mechanism. 
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Figure  17-5.  Transformation  of 
Polar  Coordinates  R  (Earth  Radii) 
and  X  (Degree)  to  Parameters  Bn, 
and  L.  (From  Mcllwain,  1961) 
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Figure  17-6.  Maximum  Proton 
Energy  aa  a  Function  of  the  Mag¬ 
netic  Shell  Parameter,  Calculated 
(Tom  Breakdown  of  the  Trapping 
Mechanism  (upper  curve)  and  from 
Explorer  XI  Observe  tiona  (lower 
curve).  (From  Germ  ire,  1963) 
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In  the  outer  zone  the  loss  processes  of  low  energy  protons  and  electrons  are 
complex.  During  observations,  the  measurement  position  varies  with  time  in  such 
a  way  that  it  's  difficult  to  separate  time  variations  from  spatial  variations. 
Nevertheless,  adiabatic  and  non-adiabatic  accelerations,  diffusion,  and  pitch-angle 
scattering  have  been  identified  as  processes  which  cause  time  variations  in  the 
outer  region. 

17.2.3  Effects  of  Nuclear  Weapons  Tests 

The  detonation  of  nuclear  weapons  at  high  altitude  injects  electrons  into  the 
trapping  region  and  also  causes  a  re-distribution  of  the  trapped  particles  that  were 
present  in  the  natural  environment.  The  largest  perturbations  observed  were  pro¬ 
duced  by  the  megaton  weapons  tested  in  1962. 

Table  17-1  gives  some  of  the  important  parameters  associated  with  some  ol 
these  nuclear  tests.  The  apparent  mean  lifetime  (the  time  to  decay  to  l/e  of  the 
initial  intensity)  of  the  injected  electrons  differs  radically  from  test  to  test.  The 
Starfish  burst  (USA  test  on  9  July  1962)  injected  electrons  with  a  long  mean  life¬ 
time;  these  electron  -  were  injected  into  a  region  of  stable  trapping  (L  **  1.  2),  and 
some  were  injected  in  equatorial  regions  at  pitch  angles  near  90°.  Electrons  from 
USSR  tests  had  relatively  short  lifetimes  (they  were  not  detectable  after  6  months); 
these  tests  injected  electrons  at  higher  latitude  and  correspondingly  larger  L  values. 
The  number  and  spatial  distribution  of  the  naturally  trapped  electrons  was  not  well 
known  prior  to  the  nuclear  tests,  hence  it  is  difficult  to  firmly  establish  what 
fraciion  of  the  electrons  subsequently  observed  were  actually  injected  by  the  weapons 
tests. 

It  is  assumed  that  proton  fluxes  measured  after  nuclear  weapons  tests  are  en¬ 
tirely  natural,  because  no  source  of  high-energy  protons  is  expected  from  a  ther¬ 
monuclear  detonation.  However,  the  spatial  distribution  of  high-energy  trapped 
protons  was  affected;  the  magnitude  of  this  redistribution  is  as  yet  unknown. 
Measurements  by  Filz  and  J.oleman  (1965)  showed  a  definite  increase  in  55-MeV 
proton  flux  at  low  altitudes  following  Starfish.  These  observations  were  consistent 
with  a  pitch-angle  redistribution  of  3.  5  degrees  for  55-MeV  protons  in  the  inner 
zone.  The  larger  increases  occurred  at  the  lower  altitudes.  This  redistribution 
could  have  been  caused  by  a  magnetic  disturbance  from  Starfish,  as  sufficient 
energy  was  available,  but  the  actual  mechanism  is  not  understood.  Other  strange 
features  are  present  in  the  measurements  made  after  Starfish.  The  secondary 
peak  of  high-energy  protons  at  the  equator  (see  Figs.  17-10  and  17-14)  is  dis 
appearing,  but  it  will  be  necessary  to  make  careful  measurements  for  many  years 
to  separate  the  effects  due  to  nuclear  weapons  tests  from  the  effects  due  to  geo¬ 
magnetic  storms  a  d  solar-cycle  variations. 
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Tanle  17-1.  Electrons  Injected  into  Trapping  Zone  by  High-Altitude 
Nuclear  Tests,  1958-1962.  {From  Van  Allen,  1966) 


Burst 

Nominal 

Altitude 

Yield 

Designation 

Date 

B(km) 

(TNT  Equivalent) 

Teak 

1  Aug.  1958 

-  75 

1C  megaton 

Orange 

12  Aug.  1958 

—  45 

10  megaton 

Argus  I 

27  Aug.  1958 

—  200 

1.4  kiloton 

Argus  II 

30  Aug.  1958 

-250 

1.  1  kiloton 

Argus  III 

6  Sept.  1958 

*.480 

1. 4  kiloton 

Starfish 

9  July  1962 

-400 

1. 4  megaton 

U.S.S.R,  I 

22  Oct.  1962 

? 

Submegaton 

U.S.S.R.  II 

28  Oct.  1962 

? 

Submegaton 

u.s.s.r.  m 

1  Nov.  1962 

? 

Megaton 

Electron 

L-value 

Apparent 

Flux- 

of 

Mean 

Designation 

i  -2-1. 

(no.  cm  s  ) 

Burst 

Lifetime 

Teak 

10? 

1.  i 

Few  days 

Orange 

103 

1.  1 

Few  days 

Argus  I 

103 

1.7 

3  weeks 

Argus  II 

105 

2.  1 

3  weeks 

Argus  III 

106 

2.0 

1  month 

Starfish 

109 

1.  12 

1.  5  years 

U.S.S.R.  I 

107 

1.9 

1  month 

U.S.S.R.  II 

107 

2.0 

1  month 

U.S.S.R.  Ill 

107 

1.8 

1  month 

Maximum  omnidirectional  intensity  at  time  of  injection. 


The  Starfish  electrons  injected  at  high  altitude  and  low  L-value  decayed  as 
predicted  from  a  diffusion  theory  developed  by  Walt  (1966).  Figure  17-7  shows 
the  theoretical  prediction  and  experimental  data.  At  high  values  of  B  (low  altitudes) 
the  agreement  is  poor.  This  is  because  the  measurements  include  the  increased 
proton  flux  caused  by  Starfish.  (The  Geiger  counter  used  to  obtain  the  counting 
rates  did  not  distinguish  electrons  from  protons. )  Figure  17-8  shows  a  comparison 
between  experimental  and  calculated  decay  times  of  artificially  injected  electrons. 
Figure  17-9  shows  the  increase  in  electrons  of  energy  greater  than  1  MeV  injected 
by  the  USSR  tests  in  1962.  For  further  details,  see  West  (1966)  and  Katz  et  al 
(1964). 
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Figure  17-7.  Predicted  values  (solid 
curve)  for  Decay  of  Trapped  Electrons 
Injected  by  the  Starfish  Burst  Compared 
with  Observed  Counting  Rates;  the  Lower 
Value  of  a  Set  of  Experimental  Points  is 
the  Counting  Rate  Corrected  for  the  En¬ 
hanced  Proton  Background.  (From  Walt 
and  Newkirk,  1966) 


Figure  17-8.  Apparent  Mean  Lifetime 
(t)  of  Trapped  Electrons;  Decay  was  As 
sumed  to  be  of  the  Form  exp(-time/T). 
(After  Van  Allen,  1966) 
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Figure  17-9.  Change  in  Omnidirectional 
Flux  of  Electrons  of  Energy  >  1  MeV 
Observed  Before  and  After  USSR  Nuclear 
Tests  in  October  1962;  B0  is  the  Magni¬ 
tude  of  the  Magnetic  Field  at  the  Equatorial 
Crossing  Point  of  the  Particles 
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17.2.4  Trapped  Protons 

Protons  with  energies  from  thermal  to  about  one  GeV  can  be  trapped  in  the 
earth's  magnetic  field.  In  general,  the  maximum  energy  of  a  proton  trapped  in  a 
specific  region  of  the  magnetosphere  is  a  function  of  the  minimum  field  strength 
encountered  along  the  particle's  trajectory;  the  greater  the  value  of  the  minimum 
field  encountered,  the  greater  the  maximum  energy  of  the  protons.  Consequently, 
the  highest  energy  trapped-protons  occur  near  the  earth  and  are  a  principal  com¬ 
ponent  oi  the  inner  zone.  In  the  outer  zone  the  protons  have  lower  energies  be¬ 
cause  the  L-shells  on  which  they  travel  may  extend  to  many  earth  radii,  and  the 
minimum  value  of  the  geomagnetic  field  strength  is  correspondingly  small.  The 
peak  flux  of  the  high  energy  protons  is  nearer  the  earth's  atmosphere,  while  peaks 
for  lower  energies  are  progressively  farther  out.  Figures  17-10  and  17-11 
illustrate  this  variation  outward  along  the  equator. 

For  protons  with  energies  above  10  MeV,  the  flux  is  relatively  stable  with 
time.  Good  model  environments  have  been  compiled  by  Vette  (1966),  which  are 
adequate  for  most  purposes.  These  proton  models,  however,  do  not  describe  the 
temporal  variations  and  are  not  completely  accurate  at  low  altitudes. 

The  processes  affecting  the  behavior  of  trapped  protons  are  the  same  as 
those  affecting  the  trapped  electrons  (Sec.  17.  2.  5).  In  the  inner  zone  (excepting 
the  lower  edge  where  particle  lifetime  is  controlled  by  the  atmospheric  density), 
proton  fluxes  tend  to  be  stable  in  time  and  fairly  insensitive  to  magnetic  activity. 

In  the  outer  zone  beyond  L«2,  decreases  in  the  proton  flux  have  been  associated 
with  magnetic  storms.  Mcllwain  (1965)  describes  a  sharp  decrease  in  proton  flux 
which  occurred  in  less  than  a  day.  This  decrease  was  energy  dependent  on  a  given 
L-shell,  but  the  mechanism  for  this  type  of  variation  is  not  yet  understood. 

Figure  17-12  gives  proton  energy  spectra  measured  at  various  L-shells. 
Figure  17-13  shows  similar  spectra  for  locally  mirroring  protons  at  the  magnetic 
equator.  These  data  indicate  there  is  no  typical  energy  spectrum;  the  form  of 
the  spectrum  varies  drastically  with  location.  Figures  17-10  and  17-14  show  the 
secondary  peak  in  high-energy  proton  fluxes;  it  is  not  known  at  this  time  whether 
this  secondary  peak  is  a  natural  feature  or  a  perturbation  resulting  from  high 
altitude  nuclear  weapons  tests.  Detailed  comparison  of  counting  rates  from  Relay  1 
shows  that  the  softest  spectrum  occurs  at  L«l.  9  a*..  the  spectra  become 
harder  toward  higher  L  values  (Fillius  and  Mcllwain,  1964). 

Figure  17-15  gives  the  distribution  of  high-energy  protons  (60  to  120  MeV) 
in  geographic  longitude  and  latitude  over  the  South  Atlantic  at  about  670-km  altitude. 
Figure  17-16  shows  the  energy  spectrum  of  trapped  protons  at  low  altitude 
measured  before  and  after  the  Starfish  nuclear  test.  The  agreement  indicates  that 
the  shape  of  the  energy  spectrum  did  not  change  drastically. 
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Figure  17-10.  Directional  Intensities 
of  Protons  Mirroring  at  the  Mag¬ 
netic  Equator;  Observations  from 
Relay  I  (Dec.  1962-March  1963) 
Compared  to  Explorer  XV  (Nov. 

1962 -Feb.  1963).  (After  Fillius 
and  Mcllwain,  1964) 


Figure  17-11.  Directional  Intensities  of 
Protons  with  Equatorial  Pitch  Angle  60°. 
(After  Davis  and  Williamson,  1966) 
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Figure  H-H.  Isointenaity  Contour*  for  40  to  110  MeV  Proton* 
(Explorer  XV  Observation*.  Nov.  1962-Feb.  1963).  Plotted  in 
Polar  Coordinate*  R  and  X;  Number*  Given  on  Contour*  are  Unit* 
of  Login  (1.4  -o).  where  Jo  is  the  Omnidirectional  Flux  (Proton* 
cm*2*~l).  The  secondary  peak  at  L  "*  2. 2  wa*  found  to  be  a 
transitory  feature;  it  may  be  an  effect  of  the  1962  nuclear  weapon*. 
(Alter  Mcllwain,  1963) 
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Figure  17-15.  Isointensity  contours  for  Protons  Trapped  Over  the  South  Atlantic 
at  Altitude  670  ±  15  km.  Fluxes  given  on  contours  are  protons  cm'2s'  *;  the 
smaller  numbers  are  typical  altitudes  (km)  at  which  measurements  were  made. 
Dashed  lines  are  contours  of  constant  L-value.  (From  Freden  and  Paulikas.  1964) 
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Figure  17-16.  Omnidirectional  Energy 
Spectrum  of  Protons  at  406-km  Altitude 
Over  the  South  Atlantic.  Circles  indicate 
data  taken  several  months  after  Starfish, 
at  L  •  1.4  and  B  •  0.218  gauss.  Squares 
Indicate  data  prior  to  Starfish;  filled 
squares  for  L  •  1. 47,  open  squares  for 
L  *  i.  35.  (From  Heckman  and  Nakano, 
1965) 
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Figure  17-17.  Intensities  of  Locally  Mirroring  Protons  In  Various  Magnetic 
Shells;  Dashed  Lines  are  Extrapolated  or  Interpolated.  The  magnetic  equator  is 
the  minimum  value  of  for  the  particular  L  value  except  for  L  ■  2.  1  and  1.  9 
earth  radii;  these  are  not  extrapolated  to  their  minimum  Bm»  (From  Fillius  and 
Mcllwain.  1964) 


Figure  17-17  shows  observations  of  the  directional  Intensity  of  low-energy 
protons;  the  flux  of  locally  mirroring  protons  is  plotted  as  a  function  of  the  mag¬ 
netic  field  at  the  mirror  point  for  various  magnetic  shells.  The  intensity  and  the 
spatial  distribution  of  these  low-energy  protons  are  not  in  accord  with  predictions 
based  on  the  hypothesis  of  a  neutron-albedo  source. 

The  most  extensive  time  histories  of  energetic  protons  in  the  lower  edge 
of  die  inner  sone  are  available  from  recoverable-satellite  nuclear-emulsion 
measurements  (Fils  and  Holeman,  1966;  Fils.  1966).  Figure  17-18  shows  the 
S8-MeV  flux  measurements  from  1961  to  1967  normalised  to  altitudes  of  278, 

350  and  400  tun.  The  solid  curves  following  July  1962  show  the  theoretical  decay 
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I  Igurc  17*18.  rime  Histories  of  5S-Mev  Fro  ion  Fluxes  at  Varioua  Attitudes.  Fhc 
solid  curves  represent  the  theoretical  decay  expected  following  Surfish  (July 
p  (atomic  electrons  per  cubic  centimeter)  is  the  measure  of  average  atmospheric 
density  used  to  compute  the  theoretical  decay,  rhe  dashed  lines  show  the  increase 
caused  by  Starfish,  the  da«h*dot  curves  show  the  time  behavior  predicted  on  the 
basis  of  solar*cycle  changes  in  atmospheric  density  and  a  neutron -albedo  source 
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based  on  ionization  loss  in  the  average  atmosphere  (Cornwall  et  al,  1965).  is 
the  average  number  of  atomic  electrons  per  cubic  centimeter  of  atmosphere  at 
the  given  altitudes,  and  is  used  as  a  measure  of  the  mean  atmospheric  density  en¬ 
countered  by  trapped  protons.  The  predicted  time  variation  of  proton  fluxes  due  to 
solar-cycle  changes  in  atmospheric  density  (Blanchard  and  Hess,  1964)  is  shown 
at  330  and  -150  km  in  F’igure  17-18.  Figure  17-19  gives  the  altitude  distribution  of 
55-MeV  protons  observed  prior  to  Starfish.  Figure  17-20  compares  these 
measurements  with  the  distributions  predicted  by  the  neutron -albedo  hypothesis 
of  particle  injection;  the  predicted  values  are  multiplied  by  a  factor  of  50  to  obtain 
the  plotted  curves.  The  measurements  of  55-MeV  protons  are  representative  of 
high-energy  protons  in  the  inner  radiation  zone.  The  increase  in  proton  flux  fol¬ 
lowing  Starfish  has  prevented  precise  measurements  of  the  natural  variation  with 
time,  of  proton  fluxes  in  the  lower  region  during  the  period  of  observation. 

The  measurements  do,  however,  lead  to  two  important  conclusions.  First, 
the  decay  of  protons  after  Starfish  is  consistent  with  that  expected  from  energy 
loss  by  ionization  in  the  upper  atmosphere.  Second,  the  decay  levels  off  as  ex¬ 
pected,  because  at  these  altitudes  the  atmospheric  density  was  decreasing  as  the 
11-year  sunspot  cycle  approached  minimum  in  October  1964.  When  corrected  lor 
the  Starfish  residual,  the  1965  fluxes  are  about  twice  the  1962  fluxes  (Filz,  1967). 
This  is  reasonably  consistent  with  the  calculations  of  Blanchard  and  Hess  (1964). 
Calculations  of  fluxes  using  the  neutron-albedo  hypothesis  do  not  agree  with  the 
observations.  As  Figure  17-20  shows,  the  undisturbed  fluxes  measured  in  1961- 
1962  are  larger  by  a  factor  of  fifty  than  those  predicted  by  calculations  of  the 
neutron-albedo  source  strength. 

1 7.2.5  Trapped  Electrons 

The  temporal  behavior  and  the  spatial  extent  of  trapped  electrons  with 
energies  above  40  keV  depend  upon  geomagnetic  activity  and  upon  solar  activity 
that  influences  conditions  in  the  magnetosphere.  Geomagnetic  activity  results  in 
(and  can  be  the  result  of)  significant  changes  in  the  flux  and  spectra  of  the  trapped 
electrons.  Figure  17-21  shows  the  correlation  between  the  variation  in  geomagnetic 
activity  (as  denoted  by  the  magnetic  activity  index  Kp)  and  the  variations  In  flux  on 
different  L-shells  in  the  outer  zone.  Figure  17-22  illustrates  the  variation  in 
spectra  with  geomagnetic  activity.  Because  the  spectra  obtained  in  the  outer  region 
are  highly  variable,  there  is  no  typical  electron  spectrum.  Sometimes  the  spectra 
depend  upon  the  L-value,  as  expected  for  particles  whose  source  is  radial  diffusion 
conserving  the  first  two  adiabatic  invariants.  At  other  times  the  spectra  are  in¬ 
dependent  of  the  L-value.  The  model  electron  environment  developed  by  Vette 
et  al  (1966)  is  the  best  model  available  at  this  time;  this  model  is  for  inner  and 
outer  zone  electrons  with  energies  greater  than  0.  5  MeV. 
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Figure  17-19.  Altitude  Distribution  of 
55-MeV  Protons  During  the  Period 
August  1961  to  June  1962  in  the  Region 
of  the  South  Atlantic  at  35°S  Latitude. 
The  altitudes  given  are  altitudes  of  the 
guiding  center;  the  solid  line  is  the 
weighted  least-squares  fit  to  the  data, 
(From  Filz  and  Holeman,  1965) 


Figure  17-20.  Omnidirectional  Fluxes  of  55-MeV 
Protons  Determined  by  Experiment  Compared 
with  Those  Predicted  by  Calculations  of  Neutron- 
Albedo  Source  Strengths;  Calculated  Values  were 
Multiplied  by  50  to  Obtain  the  Dashed  Curves. 
(After  Dragt  et  al,  1966) 
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Figure  17  -2 1.  Decay  of  Electrons 
with  Energies  Greater  than  0. 5  MeV 
at  the  Geomagnetic  Equator  (B  *  B0) 
in  the  Outer  Zone  (L  Shells  4. 0  Earth 
Radii  and  5. 0  Earth  Radii)  Compared 
with  the  Variation  in  Magnetic  Ac¬ 
tivity  Index,  Kp.  (.After  C.  E. 
Mcllwain,  "Processes  Acting 
upon  Outer  Zone  Electrons:  1, 
Adiabatic  Perturbations",  Report 
USCD-SF-66-5.  1965) 


Figure  17-22.  Spectra  of 
Electrons  Trapped  in  Outer 
Zone  at  L  *  4.  0  Earth  Radii 
for  a  Geomagnetically  Dis¬ 
turbed  Day  (29  Sept.  1964), 
a  Somewhat  Disturbed  Day 
(20  Oct.  1964)  and  a  Quiet 
Day  (20  Sept.  1964);  the 
Dashed  Line  indicates  an 
Upper  Limit  to  the  flux. 
(From  Pfitzer  et  al,  1966) 


Figure  17-23.  Typical  Energy  Spectra  of  Electrons 
in  Various  L  Shells  (1.  3  to  3.  0  Earth  Radii);  Flux 
Values  Include  Electrons  Injected  by  the  1962  Nuclear 
Weapon  Tests.  (From  Pfitzer  et  al,  1966) 
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Figure  17-24.  Contours  of  Constant  Intensity  of  Inner  Zone  Electrons 
with  Energies  Greater  than  1. 0  MeV  ±  0.  1  MeV  During  July -August 
1963;  Flux  Values  Given  are  Electrons  cm~28-lsr_1.  (Data  from 
Paolini  and  Theodoridis,  1967) 

The  inner  zone  is  much  more  stable  than  the  outer  zone,  hence  typical 
energy  spectra  and  intensity  maps  are  available.  Figure  17-23  gives  the  omni¬ 
directional  differential  spectra  for  electrons  during  November  1962  at  various  L- 
values  in  the  inner  zone.  Because  these  measurements  include  electrons  injected 
by  the  1962  nuclear  weapons  tests  and  fission-produced  electrons  have  a  character 
istic  hard  spectrum,  the  spectra  in  Figure  17-23  should  be  considered  an  upper 
limit.  Figure  17-24  is  a  typical  contour  map  for  electrons  with  energies  greater 
than  1  MeV  obtained  during  August  1963. 

The  processes  responsible  for  time  variations  in  electrons  trapped  in  the 
outer  region  may  be  grouped  phenomenologically  into  four  categories:  (1)  rapid 
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Figure  17-25.  Apparent  Rate  of 
Inward  Radial  Diffusion  of  Electrons 
with  Energies  Greater  than  1.  6  MeV; 
the  Straight  Line  has  been  Drawn 
Through  the  Data  Points .  (After 
Frank.  1966  and  Brown.  1966). 


non-adiabatic  accelerations  resulting  in  the  gain  (or  loss)  of  particle  number  or  of 
particle  energy;  (2)  adiabatic  accelerations  with  a  gain  (or  loss)  of  particle  energy; 
(3)  non-adiabatic  diffusion;  and  (4)  persistent  decay.  At  any  particular  time  one 
process  may  dominate,  although  all  processes  can  occur  simultaneously. 

Rapid  non-adiabatic  acceleration  of  electrons  may  be  caused  by  magnetic  and 
electric  fields  which  depend  on  time  and  longitude;  such  fields  can  result  from 
plasma  instabilities  (Mcllwain,  1966;  Chang,  1966;  and  Kennel  and  Petschek,  1966). 
When  non-adiabatic  acceleration  occurs,  there  may  be  changes  in  the  number  and/or 
the  energy  of  the  electrons.  The  sharp  increases  in  the  omnidirectional  fluxes 
shown  in  Figure  17-21  are  probably  due  to  a  non-adiabatic  acceleration. 

The  magnetic  fields  due  to  ring  currents  and  currents  on  the  boundary  of  the 
magnetosphere  may  result  in  adiabatic  accelerations  of  outer-zone  electrons.  The 
changes  in  electron  intensity  caused  by  ring  currents  can  be  identified  in  the  ex¬ 
perimental  data  (Mcllwain,  1966).  The  time  variations  of  these  magnetic  fields 
result  in  betatron  acceleration  of  electrons  which  are  already  trapped. 
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Figure  17-26.  Electron  Energy  Spectra  in  the  Lower  Region  of  the  Inner 
Belt  Over  the  South  Atlantic  (6 1°W  Longitude,  22°  to  35°  S  Latitude)  on 
30  Oct.  1963  (1733  to  1738  hour  Universal  Time);  Data  Points  are  not 
Corrected  for  Detector  Response.  The  fluxes  shown  here  are  the  result 
of  a  selective  redistribution  and  are  higher  than  those  observed  on  prior 
satellite  passes  in  the  same  region.  (From  Imhof  and  Smith,  1966) 


Diffusion  of  electrons  across  L-shells  is  the  result  of  a  breakdown  of  the  third 
adiabatic  invariant  while  the  first  two  invariants  are  conserved  (Dungey,  1965). 
Under  these  conditions,  electrons  gain  energy  as  they  diffuse  radially  inward;  the 
diffusion  rate  will  depend  upon  the  L-value.  Figure  17-25  shows  the  apparent  rate 
of  inward  radial  diffusion  for  electrons  measured  by  Explorers  IV,  XIV,  and  XV. 
Diffusion  across  L-shells  can  also  be  caused  by  what  is  known  as  L-shell  splitting 
and  pitch-angle  diffusion  (Roederer,  1967).  Because  of  the  day-night  asymmetry 
of  the  magnetosphere,  particles  on  the  same  L-shell  mirror  at  different  points 
on  the  night  side  of  the  earth  than  on  the  day  side.  When  particles  undergo  pitch- 
angle  changes  on  the  day  side,  they  change  both  their  mirror  point  and  L-value. 
This  type  of  diffusion  does  not  significantly  change  the  particle  energy. 

Electrons  trapped  in  the  outer  zone  undergo  a  persistent  decay.  Figure  17-21 
shows  examples  of  the  decay  on  two  different  L-shells.  This  persistent  decay  of 


31 


the  omnidirectional  flux  on  the  geomagnetic  equator  is  approximately  exponential 
in  time,  as  the  semi-log  plot  of  Figure  17-21  demonstrates.  If  the  dominant  loss 
mechanism  removes  a  fixed  fraction  of  electrons  per  unit  time,  then  the  decay  of 
the  observed  flux  will  be  of  the  form  exp  (-time/v);  r  is  observed  to  be  of  the  order 
of  two  weeks.  Non-adiabatic  pitch-angle  scattering  into  the  loss  cone  by  inter¬ 
actions  with  electromagnetic  waves  and  the  subsequent  loss  of  particles  in  the 
atmosphere  is  consistent  with  this  decay  process. 

The  processes  which  modulate  electrons  trapped  in  the  outer  zone,  influence 
the  electrons  trapped  in  the  inner  zone  to  a  lesser  extent.  For  example,  selective 
redistribution  of  electrons  has  been  observed  at  low  altitudes;  this  may  result 
from  a  combination  of  the  above  phenomena.  Figure  17-26  illustrates  the  result 
of  a  selective  redistribution  on  the  electron  energy  spectra  in  the  lower  region  of 
the  inner  zone.  The  enhanced  fluxes  shown  in  Figure  17-26  were  probably 
caused  by  a  redistribution  of  electrons  in  the  inner  zone.  This  redistribution  was 
correlated  quite  strongly  with  magnetic  activity. 

The  radial  diffusion  coefficient  for  inner  belt  electrons  has  been  measured  by 
observing  artificially  injected  electrons  (Brown,  1966).  As  Figure  17-25  shows, 
the  diffusion  coefficient  is  much  smaller  in  the  lower  edge  of  the  inner  zone  than 
’s  characteristic  of  the  rest  of  the  inner  zone. 


17.3  COSMIC  RADIATION 

Primary  cosmic  radiation  is  a  small  flux  of  high-energy  particles  of  extra¬ 
terrestrial  origin.  Galactic  cosmic  radiation  is  of  galactic  or  extragalactic  origin 
and  is  believed  to  exist  throughout  all  space  unoccupied  by  dense  matter.  The  flux 
of  galactic  cosmic  radiation  is  essentially  isotropic  with  some  minor  deviations. 
The  galactic  cosmic-ray  flux,  as  observed  at  the  earth,  is  modulated  by  the  inter¬ 
planetary  magnetic  field,  and  is  minimum  at  periods  of  maximum  rolar  activity. 

The  primary  cosmic  radiation  that  reaches  the  earth  is  composed  of  galactic 
and,  at  times,  solar  particles.  Within  the  magnetosphere  these  particles  follow 
curved  paths  dictated  by  the  magnetic  field.  The  primary  particles  incident  on  the 
top  of  the  atmosphere  collide  with  nuclei  in  the  atmosphere  initiating  a  nuclear 
cascade  of  secondary  cosmic-rays;  the  ratio  of  primary  to  secondary  particles  is 
a  function  of  altitude. 

17.3.1  Conpoaltion  and  Energy 

The  galactic  cosmic-ray  flux  is  83%  protons,  12%  alpha  particles,  1%  nuclei 
of  other  elements  (often  called  heavy  primaries),  1%  gamma  rays  and  3%  electrons 
(about  1/10  of  the  electrons  are  positive).  Antimatter  has  not  yet  been  identified 
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in  the  primary  flux;  if  it  is  present,  its  flux  must  be  less  than  0.  01%  (present 
experimental  limit)  of  the  primary  flux.  Figure  17-27  shows  relative  cosmic-ray 
abundances  of  primary  cosmic-ray  nuclei  compared  with  universal  abundances  of 
the  elements.  The  much  greater  proportion  of  lithium,  beryllium,  and  boron 
present  in  the  cosmic  radiation  is  attributed  to  their  production  by  fragmentation 
of  heavier  cosmic-ray  particles  colliding  with  the  interstellar  gas. 

The  background  flux  of  primary  protons  and  alpha  particles  impinging  on  the 

top  of  the  atmosphere  with  energies  greater  than  500  MeV  is  about  0.  2  proton 

-2-1-1  -2-1-1 

cm  s  sr  and  0. 03  alpha  particle  cm  s  sr  respectively.  Figure  17-28 

presents  the  best  available  rigidity  spectra  for  protons,  alpha  particles  and  heavy 
nuclei  of  charge  greater  than  6  at  solar  minimum.  Figure  17-29  gives  the  dif¬ 
ferential  energy  spectra  for  protons  and  alpha  particles,  also  at  solar  minimum. 

Accurate  data  for  galactic  electrons  are  not  available  at  this  time.  Figure 
17-30  summarizes  the  current  state  of  knowledge;  both  experimental  and  theoretical 
data  are  combined  to  yield  the  best  differential  energy  spectra  available. 

Figure  17-31  gives  the  integral  energy  spectrum  and  shows  the  available  experi¬ 
mental  data  points. 

17.3.2  Geomagnetic  Field  Effects 

The  cosmic-ray  particles  that  arrive  at  any  point  on  the  earth's  surface  have 
undergone  deflection  in  the  geomagnetic  field;  this  has  the  effect  of  causing  ampli¬ 
tude  and  time  variations  in  the  observed  cosmic-ray  flux.  In  order  to  relate 
variations  observed  on  the  ground  to  variations  of  the  primary  cosmic  radiation 
in  space,  it  is  essential  to  make  allowance  for  these  deflections. 

In  studies  of  time  variations  of  cosmic  radiation,  a  knowledge  of  the  dependence 
of  the  detector  counting  rate  on  the  asymptotic  direction  (the  direction  from  which 
the  particle  was  traveling  before  being  deflected  by  the  geomagnetic  field)  is 
necessary.  In  such  studies  the  concept  of  the  asymptotic  cone  of  acceptance  is 
used.  This  may  be  defined  as  the  solid  angle  containing  those  asymptotic  directions 
of  approach  which  make  a  significant  contribution  to  the  counting  rate  of  the 
detector  (or  more  simply,  the  acceptance  solid  angle  in  celestial  coordinates). 
McCracken  et  al  (1965)  and  Shea  et  al  (1968a)  list  the  geographic  coordinates  and 
the  asymptotic  directions  of  approach  for  most  of  the  cosmic -ray  stations  around 
the  earth. 

Variations  of  the  cosmic -ray  intensity  with  latitude  (and  to  a  lesser  extent 
with  longitude)  are  due  to  differences  of  the  geomagnetic  cutoff  rigidity  at  various 
locations  on  the  earth.  If  the  rigidity  (see  Eq.  17-8)  of  the  primary  nucleus  is 
equal  to  or  greater  than  the  cutoff  rigidity  for  a  given  location,  the  particle  can 
penetrate  the  geomagnetic  field  and  arrive  at  the  specified  location.  The  classic 
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Figure  17-27.  Relative 
Abundance  of  Nuclei  in  the 
Primary  Galactic  Cosmic 
Radiation  (Dashed  Line) 
Normalized  to  the  Cosmic 
Abundance  of  Hydrogen,  and 
Relative  Cosmic  Abundance 
of  the  Elements  (Solid  Line). 
(After  C.J.  Waddington. 
"Progress  in  Nuclear  Physics, 
J,  G.  Wilson,  ed. .  v.  8  Per- 
gammon  Press,  N.  Y. ,  1860) 


MMlOITY  (GV> 


Figure  17-28.  Rigidity 
Spectra  of  Galactic  Cosmic 
Ray  Particles  near  Solar- 
Cycle  Minimum:  P,  Protons; 
a,  alpha  particles;  H,  Nuclei 
with  Charge  Greater  than  6. 
(After  Anand  et  al,  1868) 


Figure  17-28.  Differential 
Energy -per -nucleon  Spectra 
of  Galactic  Cosmic  Ray 
Protons  and  Alpha  Particles 
Observed  near  die  Time  of 
Solar  Minimum.  (After 
Gloeckler  and  Jokipii,  1967) 


DIFFERENTIAL  ELECTRON  Flux  (no  m'*i‘'ftr'l6oVl) 
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Figure  17-30.  Differential  Energy 
Spectrum  of  Primary  Cosmic  Ray 
Electrons.  The  solid  line  is  a  least- 
squares  fit  to  the  data  for  energies 
above  1.  5  GeV.  The  dot-dash  curve 
represents  calculated  equilibrium 
energy-spectra  of  electrons  produced 
by  cosmic-ray  interactions  in  the 
Galaxy.  (After  L’Heureux.  1967) 


Figure  17-31.  Integral  Spectrum  of 
Primary  Electrons  of  Galactic  Cosmic 
Rays;  the  Solid  Line  is  a  Least-Squares 
Fit  to  the  Data  for  Energies  Above  1. 5 
GeV.  (After  L'Heureux,  1967) 


work  of  Stormer  on  the  motions  of  particles  in  a  dipole  magnetic  field  results  in 
the  following  equation  foi  -termining  the  cutoff  rigidity.  Pc: 

Pc  *  Me  cos4X/R^,  £l  +  (1  -  sin  e  sin  d  cos3X)*^J  (17-18) 

-5  3 

where  M  is  the  geomagnetic  dipole  moment  (3.  1165  x  10  T-  Rg  ; 

198.77  V*s*Rg  -m"1;  0.31165  gauss  *Rgj  8.05 xlO25 gauss- cm^l.c  is  the 
speed  of  light.  X  is  the  geomagnetic  latitude,  Rg  is  the  earth's  radius,  e  is  the 
zenith  angle  and  d  is  the  azimuthal  angle  measured  from  the  north.  Equation 
17-18  is  valid  in  SI  (MKSA)  and  EMU  systems  of  units. 

By  common  usage,  the  cutoff  at  a  specific  location  is  the  minimum  rigidity 
for  which  cosmic  rays  can  arrive  at  the  location  from  vertical  incidence  (zenith 
angle  zero).  For  particles  arriving  vertically,  Eq.  17  - 18  reduced  to 

Pc  *  Me  cos4X/4  R2  (17-19) 

or, 

Pc  •  14.9  cos4X  [GVJ  .  (17-20) 

The  approximation  does  not  allow  for  deviations  of  the  earth's  magnetic  field  from 
the  dipole  model  due  to  the  displacement  of  the  geomagnetic  center  from  the  geo¬ 
center,  for  deviations  due  to  magnetic  anomalies,  or  for  the  actual  magneto  spheric 
configuration.  In  a  dipole  field  cos2  X  ■  R£/L.  so  in  terms  of  the  Mcllwain  coordinate 
L,  Eq.  (17-20)  becomes 

Pc  «  14,9  L’2  JGV|  .  (17-21) 

No  analytic  equations  exist  for  calculation  of  cutoff  rigidities  in  the  actual  geo¬ 
magnetic  field,  although  many  approximations  have  been  tried.  In  practice,  a 
unique  threshold  does  not  exist  In  mid-latitude  regions  where  there  are  alternating 
bands  of  allowed  and  forbidden  rigidities  near  the  cutoff  rigidity  (penumbral  effect). 
The  method  currently  accepted  as  most  accurate  is  a  calculation  of  the  trajectories 
of  cosmic-rays  through  a  spherical  harmonic  expansion  of  the  earth's  field 
(See.  11. 2. 2. 1)  In  order  to  determine  if  a  given  particle  rigidity  is  allowed  or  for¬ 
bidden  at  a  specified  location  "^ble  17-2  lists  by  geographic  longitude  and 
latitude  the  best  available  vertical  cutoff  rigidities  determined  by  the  trajectory 
tracing  method  (Shea  et  ai,  1969b). 


Table  17-2.  Vertical  Cutoff  Rigidities  (GV) 
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Table  17-2.  Vertical  Cutoff  Rigidities  (GV)  (Contd.  ) 
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1 7.3.3  ('.harged-Particle  Flux  ftilhin  the  \imosphere 

The  intensity  and  composition  of  the  cosmic  rays  observed  within  the  atmos¬ 
phere  depend,  in  addition  to  the  cutoff  rigidity  of  the  observation  point,  also  on  the 
quantity  of  absorbing  material  traversed  before  observation.  The  ionization  rate 
measured  within  the  atmosphere  depends  upon  the  amount  of  matter  above  the 
point  of  observation  and  on  its  distribution  in  height.  Atmospheric  conditions, 
especially  barometric  pressure,  have  an  appreciable  effect  on  the  measured  in¬ 
tensity,  hence  cosmic-ray  intensities  are  usually  reported  in  terms  of  atmospheric 
depth  (mass  of  air  per  unit  area  above  the  observation  point)  or  of  air  pressure  at 
the  observation  point  rather  than  the  altitude  of  the  observation. 

At  a  given  altitude,  the  value  of  the  pressure  in  millibars  is  about  2%  less 
than  the  atmospheric  depth  in  grams  per  square  centimeter.  Figure  17-32  is  a 
plot  of  atmospheric  depth  as  a  function  of  altitude;  for  graphs  and  tables  for  com¬ 
puting  air  masses  at  various  zenith  angles,  see  Ely  (1S62).  Figures  17-33  and 
17-34  illustrate  the  variations  of  cosmic  ray  intensities  with  atmospheric  depth 
and  with  solar  activity;  ionizatic  rates  in  the  atmosphere  are  given  as  a  function 
of  atmospheric  depth  for  a  ter  -•  ar  period  that  includes  the  solar-cycle  minimum 
in  1954  and  the  maximum  in  1957.  Figure  17-35  shows  the  different  ionization 
rates  in  the  atmosphere  at  roughly  30-km  altitude  as  a  function  of  north  geomagnetic 
latitude  during  four  different  years. 

Primary  cosmic  rays  incident  upon  the  top  of  the  earth's  atmosphere  interact 
with  air  nuclei  producing  high-energy  secondary  cosmic-rays.  These  secondary 
particles,  in  turn,  interact  with  other  nuclei  and  produce  additional  secondaries. 
Figure  17-36  illustrates  a  nuclear  cascade  process  initiated  by  a  primary  cosmic 
ray.  The  production  of  secondary  components  becomes  significant  at  about  55-km 
altitude  (4-mb  pressure)  with  the  local  intensity  reaching  a  maximum  (the  Pfotzer 
maximum)  at  approximately  20  km  (56  mb).  The  intensity  of  secondaries  then 
decreases  from  the  Pfotzer  maximum  to  the  surface  of  the  earth  as  the  particles 
lose  energy  by  additional  collisions  until  the  majority  either  decay  cr  are  absorbed. 

The  particles  that  comprise  the  secondary  component  are  primarily  v  and 
n  mesons,  neutrons,  protons,  electrons,  and  y-ray  photons;  neutrons  and  /u  mesons 
are  the  secondaries  most  commonly  measured  on  the  earth's  surface.  The 
secondary  cosmic  rays  are  often  classified  into  three  major  components:  the  hard 
component,  which  is  primarily  relativistic  mesons;  the  nucleonic  component, 
which  is  the  locally  produced  protons  and  neutrons;  and  the  soft  component,  which 
is  the  electrons  and  the  y  rays.  Figure  17-37  shows  the  intensity  at  solar  mini¬ 
mum  of  fast  secondary  neutrons  at  four  north-geomagnetic  latitudes  as  a  function 
of  atmospheric  pressure. 
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Figure  17-33.  Yearly  Average  of 
Cosmic  Ray  Ionization  Rate  per 
Atmosphere  of  Air  as  a  Function 
of  Atmospheric  Depth.  (From 
Neher  and  Anderson,  1962) 
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Figure  17-34.  Cosmic  Ray 
Ionization  Rate  per  Atmos¬ 
phere  of  Air  as  a  Function  of 
Time  for  Selected  Atmospheric 
Depths.  Quarterly  Zurich 
sunspot  number  and  quarterly 
planetary  magnetic  character 
figures  are  shown  for  com¬ 
parison.  (From  Neher  and 
Anderson,  1962) 
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Figure  17-35.  Cosmic  Ray 
Ionization  Rates  per  Atmosphere 
of  Air  Near  the  Top  of  the  Atmos¬ 
phere  as  a  Function  of  Latitude  for 
Various  Years.  (From  Neher  and 
Anderson  1962) 
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Figure  17-36.  Schematic  Diagram  of  a 
Cosmic  Ray  Shower.  N  and  P  are  high 
energy,  n  and  p  disintegration-product, 
neutrons  and  protons;  pions,  mesons, 
electrons,  positrons,  and  gamma  rays  are 
indicated  by  conventional  symbols 
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Figure  17-37.  Flux  of  1  to  10  MeV  Neutrons 
as  a  Function  of  Altitude  (Atmospheric  Pressure) 
at  Various  Indicated  Latitudes  During  Solar 
Minimum.  (After  Holt,  et  al,  1966) 


In  theory,  the  intensity  of  any  specified  secondary  component  can  be  derived 
from  a  knowledge  of  the  primary  spectrum  and  the  specific  yield  functions,  but  at 
present  the  yield  functions  are  not  accurately  known.  Perhaps  the  most  advanced 
work  on  the  secondary  component  is  that  being  done  on  the  atmospheric  neutrons. 
Neutron  monitors,  which  measure  the  nucleonic  component,  are  advantageous  for 
the  study  of  primary  flux  variations  because  they  are  particularly  sensitive  in¬ 
dicators  of  primary  radiation,  with  energies  from  500  MeV  to  a  few  GeV.  The 
analysis  of  neutron  monitor  data  provides  a  great  deal  of  our  knowledge  of  cosmic- 
ray  variations;  data  for  a  complete  solar  cycle  (beginning  in  1955)  have  been  ob¬ 
tained  by  a  worldwide  network  of  neutron  monitors. 

17.3.4  Flux  Modulations 

There  is  no  experimental  evidence  that  the  average  galactic  cosmic-ray 

intensity  has  undergone  any  significant  long-period  changes.  The  radioactive 
14  3 

isotope  C  (5  x  10  yr  half-life)  is  produced  in  the  atmosphere  by  cosmic  rays, 

14 

and  agreements  in  die  dating  of  archaeological  artifacts  by  the  C  method  and  by 
archaeological  methods  imply  that  any  long-period  change  in  the  cosmic-ray 
intensity  in  the  last  30, 000  years  is  smaller  than  a  few  percent.  The  abundance 
of  Be  (2.  5  x  10  yr  half-life),  another  cosmic-ray  produced  isotope,  indicates 
that  the  mean  cosmic-ray  intensity  has  changed  by  less  than  a  factor  of  10  in  the 

C 

past  10  years.  There  are,  however,  in  addition  to  the  variation  of  intensity  during 
the  solar  cycle,  short-term  variations  in  coBmic-ray  intensities. 

The  variation  of  cosmic-ray  flux  with  the  11-year  solar  cycle  is  indicated  by 
the  year  to  year  variations  shown  in  Figures  17-33,  17-34,  17-35,  and  17-38. 
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Figure  17-38.  Monthly  Averages  of 
Relative  Intensities  Observed  by  the 
Neutron  Monitor  at  Mt.  Washington, 
N.  H. ;  Cutoff  Rigidity  1. 4 1  GV 
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An  increase  in  solar  activity  corresponds  to  a  decrease  in  galactic  cosmic-ray 
intensity.  From  solar  minimum  to  solar  maximum,  the  low-energy  primary  par¬ 
ticles  tend  to  be  removed  from  the  flux  in  the  vicinity  of  the  earth,  decreasing  the 
energy  density  of  these  primary  galactic  cosmic-rays  by  about  40%.  Figure  17-38 
illustrates  the  solar  cycle  modulation  of  the  cosmic-ray  nucleonic  intensity  ob¬ 
served  by  the  Mt.  Washington  neutron  monitor.  The  intensity  does  not  change 
smoothly  from  maximum  to  minimum  values.  It  drops  sharply  and  the  decrease 
is  followed  by  a  partial  recovery.  This  modulation  is  thought  to  be  caused  by  an 
increase  in  the  extent  of  the  solar-interplanetary  magnetic  field  contained  in 
plasmas  that  are  emitted  by  the  sun  in  increased  amounts  during  active  periods. 
Changes  in  cosmic-ray  intensity,  however,  lag  changes  in  sunspot  number  by 
9  to  12  months.  As  the  solar  activity  increases  from  minimum,  primary  particles 
of  increasingly  higher  energy  are  affected.  Particles  with  rigidity  higher  than  about 
15  GV,  however,  seem  to  be  relatively  unaffected. 

Short-term  variations  can  be  classified  as  follows:  27-day  variations;  24-hour 
variations;  sudden  decreases;  and  sudden  increases.  The  27-day  variations  are 
related  to  processes  in  the  solar  atmosphere  which  affect  conditions  of  the  inter¬ 
planetary  medium.  The  24-hour  variations  are  normally  about  0. 25%,  implying  a 
small  diurnal  anisotropy.  This  anisotropy  was  invariant  during  the  period  from 
1957  to  1965  while  solar  activity  decreased  from  maximum  (1957)  to  minimum 
(1964).  The  normal  state  of  the  diurnal  anisotropy  is  energy  independent.  The 

■  3 

amplitude  varies  as  the  cosine  of  the  declination,  and  the  maximum  flux  (3.  8  x  10 
times  the  isotropic  flux)  is  incident  from  a  direction  90°  to  the  east  of  the  earth- 
sun  line  (McCracken  and  Rao,  1966).  There  is  also  a  semi-diurnal  anisotropy 
with  an  amplitude  of  approximately  0.  1%  and  maximum  60°  to  the  west  of  the  earth- 
sun  line  (Abies,  et  al,  1966).  For  detailed  discussions  of  these  variations  and 
the  variation  with  geomagnetic  disturbances,  see  Dorman  (1963). 

Sudden  decreases,  called  Forbush  events,  are  associated  with  interactions 
between  primary  particles  and  the  magnetic  fields  contained  in  solar  corpuscular 
streams  and  interplanetary  plasma. 
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At  the  onset  of  a  geomagnetic  stcrm  (Sec.  11.  4),  or  a  few  hours  after  the 
onset,  cosmic-ray  intensity  may  decrease  rather  sharply  for  about  a  day  to  a 
minimum  that  is  sometimes  more  than  10%  below  the  pre-storm  value.  The  in¬ 
tensity  then  increases  rather  slowly,  generally  taking  several  days  to  recover, 
during  which  time  the  amplitude  of  the  diurnal  variation  may  be  enhanced. 

A  few  decreases  in  cosmic-ray  intensity  that  are  not  associated  with  visual 
solar  flares  have  been  observed  during  magnetic  storms.  Conversely,  some 
large  geomagnetic  storms  are  not  associated  with  any  significant  changes  in  the 
cosmic-ray  intensity.  Although  during  magnetic  storms  in  which  the  earth's 
field  decreases  there  is  a  corresponding  decrease  in  cutoff  rigidity,  accurate 
determinations  of  cutoff  rigidities  are  not  available  for  geomagnetically  dis¬ 
turbed  periods. 

17.4  ENERGETIC  SOLAR  PARTICLES  (SOLAR  COSMIC  RAYS) 

Charged  particles  emitted  sporadically  by  the  sun  with  energies  from  a  few 
MeV  to  a  few  GeV  are  described  here;  the  continuous  solar  emission  of  charged 
particles  with  much  lower  energies  (the  solar  wind)  is  described  in  Chapter  18. 

The  highest  energy  particles  (above  500  MeV)  are  emitted  only  at  times  of 
intense  solar  activity,  generally  during  solar  flares  of  importance  3  or  4.  The 
low  energy  particles  (about  5  to  100  MeV)  can  be  associated  with  solar  flares  of 
importance  2  or  greater.  Most  particle  events  can  be  associated  with  a  particular 
solar  flare,  but  most  solar  flares  do  not  produce  particle  events  detectable  by 
earth-based  sensors. 

Sea  level  events  (also  called  ground  level  events)  indicate  the  arrival  of  the 

highest-energy  particles.  For  example,  neutron  monitors  located  on  the  earth's 

surface  detect  the  neutrons  produced  in  the  atmosphere  by  incoming  charged 

2 

particles  of  energy  greater  than  500  MeV;  the  1000  g/cra  of  atmosphere  above 
the  monitors  absorbs  essentially  all  particles  with  less  energy  (Sec.  17.  3.  3). 
Other  types  of  ground-based  instruments  depend  on  this  absorption  for  their 
|  operation,  utilizing  the  enhanced  ionization  in  the  ionosphere  to  indicate  the  ar¬ 

rival  of  charged  particles.  Riometers,  which  are  sensitive  detectors  of  cosmic 
radio-noise  intensity,  reveal  the  influx  into  polar  regions  of  protons  with  energies 
of  about  5  to  100  MeV;  the  enhanced  ionization  produced  by  these  protons  causes 
strong  absorption  of  cosmic  radio-noise,  i.  e.  a  polar  cap  absorption  (PCA)  event. 

The  frequency  of  occurrence  and  the  intensity  of  the  events  vary  strongly  with 
the  sunspot  cycle.  As  Figure  17-39  shows,  the  number  of  PCA  events  observed 
during  the  sunspot  minimum  was  of  the  order  of  one  per  year,  whereas  during  the 
sunspot  maximum  the  rate  of  occurrence  was  roughly  one  per  month.  Table  17  -3 

i 

1 
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*  IONOSONDE  DATA,  SOME  EVENTS  MAT  HAVE  BEEN  MISSED 

Figure  17-39.  Number  of  Solar  Particle  Events  Detected  at  Sea  Level 
and  in  the  Ionosphere  (PCA  Events),  in  Relation  to  the  Solar  Cycle. 

The  Sea-Level  Events  Indicate  the  Arrival  of  the  Highest  Energy 
(>500  MeV)  Solar  Particles;  PCA  Events  Indicate  the  Arrival  of  Lower 
Energy  (> 20  MeV)  Solar  Particles.  (After  Fich tel  and  McDonald,  1967) 


describes  the  large  events  observed  during  high  solar  activity  (1956-1961). 

Table  17-4  lists  the  most  significant  events  during  low  solar  activity  (1962-1966). 
Many  small  events  were  also  observed  during  this  period  of  solar  minimum,  but 
a  direct  comparison  with  the  frequency  of  small  events  during  the  earlier  solar 
maximum  cannot  be  made  because  the  detectors  used  during  the  minimum  period 
were  much  more  sensitive  than  those  previously  available.  The  intensities  and 
durations  of  the  solar  disturbances  vary  so  much,  and  the  interplanetary  medium 
between  the  sun  and  the  earth  is  so  variable,  that  there  is  a  little  value  in  con¬ 
structing  hypothetical  models  of  an  "average"  solar  cosmic-ray  event.  This  dis¬ 
cussion  is  limited  to  describing  the  parameters  used  to  characterize  events  and 
to  indicating  the  observed  limits  of  these  parameters. 

An  energetic  solar-particle  event  may  be  characterized  by  the  following  para¬ 
meters: 

(1)  particle  types  present; 

(2)  intensity  (generally  the  number  of  particles  per  unit  area  per  second  with 
energy  above  a  specified  value); 

(3)  energy  spectrum  or  rigidity  spectrum; 

(4)  rise  time  and  decay  time; 

(5)  degree  of  anisotropy  and  its  direction;  and 

(6)  volume  of  interplanetary  space  over  which  the  flux  extends. 


Tabie  1,-3.  Some  Major  Solar  Cosmic  Ray  Observations  During  1956 
1961.  (After  McDonald  1963). 


Solar  Flare  Data 

Date 

Im¬ 

por¬ 

tance 

Helio- 

graphic 

Position 

1  deg ) 

Optical 

Maxi¬ 

mum 

fllTl 

1956  Feb  23 

3+ 

N22  W74 

0340 

1957  Jan  20 

3+ 

S25  W30 

1120 

Mar  23 

3+ 

SO  E78 

1005 

J'll  7 

3+ 

N25  W08 

0115 

1958  Aug  16 

3+ 

S14  W50 

0440 

Aug  22 

3 

N18  W10 

1448 

Aug  26 

3 

N20  WS4 

0027 

May  10 

3+ 

N23  E47 

2118 

1959  JullO 

3-f 

N22  E70 

0222 

Jul  14 

3+ 

N16  E07 

0349 

Jul  16 

3+ 

N08  W26 

2145 

Apr  1 

3 

M2  W10 

0859 

Apr  5 

2+ 

N10  W61 

0245 

Apr  28 

3 

SOS  E34 

0130 

May  4 

3+ 

N14  W90 

1020 

1960  Sep  3 

3 

N18  E88 

0110 

L  Nov  12 

11 

3+ 

N27  WQ2 

1329 

f  Nov  15 

3+ 

N30  W32 

0221 

Nov  20 

3 

N28  W113 

2020 

Jul  11 

3 

S06  E32 

1700 

Jul  12 

3+ 

SOT  E22 

1030 

i%i  jul  ia 

3+ 

S06  W60 

1010 

Jul  20 

3+ 

S07  W90 

(16001 

[  Sep  28 

C 

3 

N14  E30 

2223 

Cosmic  Ray  Data 


Onset  to  Max. 

Peak  Flux 

Integrated 

Time 

Decay  Time 

(no.  cm-2 

Intensity 

(h) 

(h) 

sec"1) 

(no.  cm-2) 

>30 

>  100 

>30 

>  100 

>30 

MeV 

MeV 

MeV 

MeV 

MeV 

6.5  X  10s  3.2  X  10* 


6.000-  1,000 


7  X  10*  7.5  X  10* 


40  20  4,000  1.200  8.8  X  10s  1.0  X  10* 


9-12  10,000-  1,200  1.1  X  10*  6.3  X  10* 


8  16,000-  1,500  8.1  X  10*  1.3  X  10* 


2.7  X  10*  15  X  10» 


2  40  2 


200  40 


32  26  240  60 


8-10  18-24  14-18  12,000  2500  1.4  X  10*  35  X  10* 


6,000  2.400  55  X  10*  1.2  X  10* 


IS  1.0  X  10*  1.6  X  10* 


600  2.1  X  10*  45  X  1C* 


25  X  10*  9.7  X  10* 


2 

4  I 

2 

2500 

•  Value*  tend  to  be  overestimated. 
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Table  17-4.  Principal  Solar  Cosmic-Ray  Events  1962-1966 


 - 

Date 

Importance 

Class 

Time  of 
Optical  Flare 
<h) 

Heliographic 

Position 

(deg) 

Peak  Flux 
Protons  >20  MeV 
(no.  cm 

1962  Feb  20* 

2 

1222 

N10  E79 

0.  5 

1963  Sep  20 

2 

2255 

N10  W09 

4 

1964  Mar  16 

2 

1550 

N05  W75 

0.7 

1965  Feb  5 

2 

1750 

N07  W25 

50 

Oct  4 

2 

0938 

S20  W29 

3 

1966  Mar  24 

3B 

0233 

N  18  W37 

15 

Jul  7 

2B 

0022 

N34  W45 

30 

Aug  28 

2B 

1530 

N08  W03 

15 

*A  number  of  flares  occurred  on  thi.^  date;  data  listed  correspond  to  the 
first  flare  occurrence. 


Unfortunately,  these  parameters  do  not  remain  constant  from  event  to  event  and 
even  change  during  a  single  event.  The  rise  and  decay  times  tend  to  be  shorter 
for  higher-energy  particles  than  for  lower-energy  ones,  so  that  generally  the 
spectra  become  steeper  at  later  times  in  an  event.  It  is  convenient  to  give  the 
time  history  of  particles  above  a  specified  threshold  energy;  because  of  the 
steepness  of  solar -particle  spectra,  this  is  essentially  the  time  history  of  the 
particles  at  the  specified  energy.  Table  17-3  gives  an  indication  of  the  time 
histories  for  30  MeV  and  100  MeV  protons  for  a  number  of  large  events. 

17.1.1  Composition 

Solar  cosmic-ray  events  are  often  referred  to  as  solar  proton  events.  Al¬ 
though  protons  are  the  major  component,  significant  numbers  of  alpha  particles 
are  observed.  The  rigidity  spectra  of  protons  and  alphas  tend  to  be  similar,  but 
the  relative  proportion  of  particles  varies  considerably  between  events  and  even 
during  a  given  event.  Within  a  given  rigidity  interval  an  alpha -to-proton  ratio  as 
high  as  1  has  been  reported,  whereas  in  some  events  very  few  alpha  particles 
have  been  detected.  This  variable  alpha-to-proton  ratio  is  one  of  the  most 
puzzling  aspects  of  solar  cosmic-rays.  Too  little  is  known  at  this  time  to  discuss 
the  abundance  of  the  various  hydrogen  and  helium  isotopes.  A  rough  estimate  for 
the  deuteron-to-proton  and  the  triton-to-proton  ratio  in  the  energy  interval  10  to 
100  MeV/nucleon  is  of  the  order  10  or  less. 


Table  17-5.  Estimated  Relative  Abundances  (By  Number  Based  on  1.  0 
for  Oxygen)  of  Solar  Particles  Compared  with  Galactic  Cosmic  Rays. 
(After  Fichtel  and  McDonald,  1967). 
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Element 

Solar 

Galactic 

Cosmic  Rays 

Cosmic  Rays 

Photosphere 

Corona 

He 

107  ±  14 

? 

445 

48 

Li 

- 

<10-^ 

- 

0.3 

Be-B 

<0.  02 

<10‘5 

- 

0.8 

C 

0.  59  ±  0.  07 

0.6 

1.3 

1.8 

N 

0.  19  ±  0.  04 

0.  1 

0.  1 

5,0.8 

O 

1.0 

1.0 

1.0 

1.0 

F 

<0.03 

0.001 

- 

5,0.  1 

Ne 

0.  13  ±  0.  02 

? 

C.  11 

0.  30 

Na 

- 

0.002 

0.  01 

0.  19 

Mg 

0.043  ±  0.011 

0.  027 

0.  20 

0.  32 

A1 

- 

0.002 

0.01 

0.  06 

Si 

0.  033  ±  0.011 

0.035 

0.22 

0.  12 

P-Sc 

0.  057  ±  0.  017 

0.032 

- 

0.  13 

Ti-Ni 

5.0.  02 

0.006 

-0.  1 

0.28 

Nuclei  with  charge  greater  than  2  have  recently  been  observed  in  a 
number  of  events.  (During  the  1954-1964  solar  cycle,  experimental  techniques 
were  not  adequate  to  establish  or  rule  out  their  presence. )  The  rigidity  spectra 
observed  tend  to  be  similar  to  those  of  the  protons  and  alpha  particles.  Table 
17-5  gives  the  best  available  data  on  the  relative  abundances.  The  relative 
abundance  of  elements  in  the  solar  cosmic  rays  is  closer  to  that  of  the  solar 
photosphere  than  to  that  of  galactic  cosmic  rays.  The  absence  of  lithium, 
beryllium,  and  boron  in  the  solar  particles  is  consistent  with  their  low  universal 
abundance.  The  admixture  of  all  elements  heavier  than  helium  1b  about  1  to  2%  of 
the  helium  component.  Apparently,  the  abundance  ratio  of  heavier  elements  to 
helium  is  reasonably  ecus  tant  as  a  function  of  particle  rigidity. 

The  el*  'tron  component  is  of  particular  interest  because  it  provides  a  possible 
explanation  for  the  observed  radio  emission  from  active  flare  regions  (Sec.  16. 3. 5). 
Calculations  (Stein  and  Ney,  1963)  satisfactorily  explain  the  radio  noise  as 
synchrotron  radiation  from  electrons  with  flux  and  initial  rigidity-spectrum 
similar  to  that  observed  for  protons.  Observed  fluxes  of  higher-energy  electrons 
are  rather  weak  although  relatively  high  inteusity  fluxes  of  low-energy  electrons 
are  observed.  A  number  of  electron  events  with  an  apparent  absence  of  protons 
have  been  detected. 


17.4.2  Isteasity  and  Speeds 

In  the  largest  events  associated  with  solar  flares,  the  flux  of  particles  near 

3  «2  - 1 

the  earth  with  energy  greater  than  30  MeV  exceeded  10  particles  cm  s  for  a 


48 


Table  17-6.  Large  Events  During  the  Period  1956  through 
1961.  (After  Fichtel  and  McDonald.  1967). 


- - - -  1 

Events 

Integrated  Flux 

E  >  30  MeV 

(no. ) 

(particles  cm*2s"l) 

10 

107  to  108 

10 

108  to  109 

2 

>  109 

period  of  more  than  a  day.  The  total  number  of  particles  arriving  at  the  earth 
during  such  events  exceeded  the  total  yearly  flux  of  galactic  cosmic  rays. 

Table  17-6  gives  the  size  distribution  of  the  largest  events  observed  from  1956 
through  1961. 


Figure  17-40.  Energy  Spectra  of  Protone  from  Several  Moderate- 
Site  Solar  Events  Compared  with  the  Galactic  Cosmic  Ray  Spectrum 
at  Solar  Minimum.  (From  Fichtel  and  McDonald.  1967) 
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MOITYttV) 


Figure  17-4 1.  Integral  Proton  Intensity  as  a  Function  of  Rigidity  at 
Selected  Times.  Data  points  taken  from  counter  ascents  are  shown 
as  solid  symbols;  those  taken  with  emulsions,  as  open  symbo's.  The 
time  of  observation  is  shown  for  each  spectrum;  the  lime  of  flares  is 
given  in  Table  17-3.  (From  Freier  and  Webber,  1963) 


Because  both  the  rise  and  decay  times  of  the  particles  vary  with  energy,  the 
energy  spectrum  is  time-dependent.  Solar  cosmic  rays  coniain  fewer  extremely 
high  energy  particles,  and  therefore  have  steeper  spectra  than  the  galactic  cosmic 
rays.  Figure  17-40  presents  a  number  of  solar  particle  spectra  in  comparison  to 
the  galactic  ..oamic-ray  spectrum.  Reasonably  successful  attempts  have  been  made 
to  express  the  time  change  in  the  differential  rigidity  spectrum  by  the  equations 
of  the  form: 

J(P)  •  Jo<0  exp  |  -P/P0(t)J  (17-22) 


PQ  It*  normally  in  the  range  40-400  MV  and  decreases  with  time  from  the  begin¬ 
ning  of  the  event.  Figure  17-4 1  shows  some  exponential  rigidity  spectra  observed 
during  various  times  within  six  events. 


Do 

17.4.3  \nisrtropy 

Anisotropy  of  solar-produced  particles  is  related  most  strongly  to  the 
storage  and  propagation  properties  of  the  magnetic  fields  in  interplanetary  space. 
Although  data  are  still  too  rudimentary  to  uniquely  determine  details  of  the 
processes  involved,  a  broad  model  has  developed. 

Analyses  of  the  results  from  the  neutron-monitor  network  show  that  early  in 
the  large  events  the  solar  particles  with  energies  greater  than  500  MeV/nucleon 
are  generally  anisotropic.  Later  on,  however,  the  flux  becomes  isotropic,  which 
suggests  a  temporary  storage  of  particles  in  space  near  the  earth.  During  the 
anisotropic  phase,  the  particles  appear  to  arrive  at  the  earth  from  a  preferential 
direction  approximately  50  degrees  west  of  the  earth-sun  line.  This  is  consistent 
with  the  concept  of  particles  flowing  away  from  the  sun  in  paths  along  solar  mag¬ 
netic  field  lines  which  have  an  Archimedes  spiral  (garden  hose)  configuration  (see 
Chapter  18).  This  is  also  consistent  with  the  observations  of  optical  flares 
associated  with  solar  cosmic-ray  events  on  the  earth;  these  flares  are  located 
more  frequently  in  the  western  hemisphere  of  the  solar  disk. 

Measurements  from  satellites  and  probes  of  solar  particles  with  energies 
from  1  to  70  MeV  and  of  interplanetary  magnetic  fields  show  that  within  the 
general  spiral  shape  the  magnetic  fields  also  have  a  finer,  filamentary  form  which 
normally  contains  some  irregularities.  These  filaments  move  with  the  solar  wind 
and  constrain  the  lower  energy  particles  to  move  along  the  filamentary  paths.  At 
these  lower  energies  the  particle  flux  can  be  highly  anisotropic  for  periods  up  to 
48  hours.  Major  changes  in  direction,  which  are  related  to  changes  in  the  mag¬ 
netic  field,  occur  from  hour  to  hour.  One  of  the  implications  of  these  observations 
is  that  there  is  some  restriction  on  the  spatial  extent  of  the  solar  particle 
propagation.  Very  little  direct  evidence  is  available,  but  it  appears  that  some¬ 
what  less  than  half  the  inner  solar  system  is  populated  by  solar  cosmic-ray  par¬ 
ticles  during  any  one  event. 
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IS.  ABSTRACT 

This  survey  is  a  summary  of  available  information  on  the  corpuscular  radiation 
environment  in  the  vicinity  of  the  earth  and  in  the  earth's  atmosphere;  radiation 
trapped  in  the  magnetosphere  (Van  Allen  belts),  galactic  cosmic  rays,  and  high- 
energy  solar  particles  (solar  cosmic  rays)  are  described.  The  topics  presented 
are:  the  motions  of  charged  particles  in  the  geomagnetic  field;  the  spatial  and 
temporal  distributions,  intensities,  and  energies  of  electrons  and  protons  trapped 
in  the  earth's  magnetic  field;  the  composition,  intensity,  energy  spectrum  and 
variations  of  galactic  cosmic  rays;  and  the  composition,  intensities,  anisotropy, 
and  variations  of  the  high-energy  charged-particles  emitted  sporadically  by  the  sun 
in  association  with  solar  flares  and  other  solar  disturuances. 
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